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I. INTRODUCTION 

Controlled manipulation of coherent quantum systems 
is a crucial requirement for quantum information tech- 
nologies, can be exploited in ultrafast switches, and may 
allow the exploration of exotic regimes of quantum dy- 
namics. An important example among solid-state sys- 
tems is that of excitons in quantum dotsP which provide 
discrete atomic-like transitions that can be manipulated 
using optical pulses. These transitions have been demon- 
strated experimentally to correspond to two-level sys- 
tems, for w hich resonant optical excitation induces Rabi 
oscillations PEl Thus under pulsed excitation the number 
of excitons created oscillates with the pulse area: a pulse 
of the correct duration and intensity creates exactly one 
exciton in the quantum dot, while other pulses create su- 
perpositions of one- and no- exciton states. Frequency- 
swept laser pulses have also been used to create single 
excitons in quantum dots, implementing the protocol of 
adiabatic rapid passage (ARP), and allowing state ma- 
nipulation in a way which is robust against fluctuations 
in the coupling strengths and transition energies of the 
dotsP3 

A theoretical description of such coherent control ex- 
periments must capture both the dynamics of the driven 
quantum dot, and the scattering and decoherence intro- 
duced by the interaction between the dot and its envi- 
ronment. In particular, the coupling to acoustic phonons 
leads to dephasin g of the Rabi oscillations^ and limits the 
inversion in ARPPEEESl (Optic phonons could play a role 
in ARP with very intense, very short pulsesp^where the 
energy scales are comparable to the optic phonon ener- 
gies, but we do not consider this regime here.) The stan- 
dard theoretical approach^ involves second-order per- 
turbation theory in the phonon coupling, leading to an 
equation of motion for the reduced density operator of 
the dot. This equation involves an integral over all pre- 
vious states of the dot, capturing the memory effects due 
to the finite bandwidth and response time of the envi- 
ronment. 



The most popular approach to treating this type of 
non-Markovian equation of motion is a form of Markov 
approximation which reduces the equation to a time- 
local equation, with the effects of the phonons appear- 
ing as a constant Lindblad form describing dephasing. 
Such an equation would be valid under the assumption 
that the response time of the environment is the short- 
est timescale in the problem. This approximation, how- 
ever, is generally incorrect for quantum dot excitons! 13 ^ 14 ! 
It predicts that the environment induces transitions in- 
dependently of the state of the dot, in contradiction to 
experiments in which the driving field changes the de- 
phasing. Such dynamical excitation-dependent dephas- 
ing has, however, been successfully described by more 
sophisticated approaches. These i nclude nu merical solu- 
tions of the full integral equation) 9 ! 10 ! 13 * 1 ^ including all 
the memory effects of th e environment, as well as time- 
local approximations^ 4 *^ to this equation, allowing for 
some of the memory effects neglected in the simplest 
form of Born-Markov approximation. This same point 
— that it is important to calculate the decay rates aris- 
ing from system-bath coupling making use of the actual 
system Hamiltonian — is also crucial in describing the 
correct equilibrium state of strongly coupled systems, as 
discussed in Ref. HH1 

In this paper, we discu ss the application of the Bloch- 
Redfield-Wangsness 19 20 (BRW) theory, widely used to 
describe nuclear spin relaxation, to coherent control ex- 
periments on excitons in quantum dots. It is similar to 
the simplest Born-Markov approximation, but allows the 
timescale set by the inverse level spacing of the Hamilto- 
nian to be smaller than the response time of the environ- 
ment. We use this theory to derive a generalized Lind- 
blad form for the phonon-induced damping, in which the 
transition operators connect the time-dependent eigen- 
states of the dot, with the expected perturbative transi- 
tion rates [Eq. Q]. This differs from the simplest Lind- 
blad form mentioned above, which is often applied to 
quantum dots, in which the transition operators do not 
necessarily connect eigenstates. Furthermore, in con- 
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trast to the fo rms o btained by some other time-local 
approximations G^llfil guarantees the positivity of the 
density operator, and so can be used across a wider va- 
riety of pulse shapes and temperatures. We focus on the 
application to determining the effects of phonons for adi- 
abatic rapid passage in quantum dots, and show that the 
results are similar to those recently obtained by numeri- 
cal solutions of the full integral equationPi^l Thus Eq. Q 
provides a simple picture of the effects of dephasing and 
a lightweight computational approach for modeling de- 
phasing in quantum dots. As a further application of 
this theory, we use it to demonstrate the feasibility of 
generalizing the ARP experiments to measure phonon 
spectra and distribution functions, implementing a type 
of multi-dimensional spectroscopy for excitons. 

The remainder of this paper is structured as follows. 
In Sec. |TlJ we outline the derivation of the equation of 
motion for a driven quantum dot interacting with acous- 



Hamiltonian may be expressed, using pseudospin-1/2 op- 
erators s, as (h — 1) 



tic phonons. In Sec. Ill we present the predictions of this 
equation for the inversion (exciton occupation) obtained 
in ARP, discuss how this process may be optimized in the 
presence of phonon-induced dephasing, and explain why 
the dephasing can in some circumstances improve the fi- 
nal inversion. In Sec. |IV| we use it to propose and analyze 
two forms of spectroscopy, in which the final inversion 
reached is used to measure the density-of-states of the 
phonons. In Sec. |Vj we discuss further the relationship 
between the positivity-preserving Lindblad form [Eq. Q] 
obtained here and the generally positivity-violating pre- 
cursor to this form [Eq. Q], which is sometimes used 
directlyP^Efil We present numerical results showing the 
failure of this latter approximation in the case of ARP 
pulses. In Sec. [Vljwe summarize our conclusions. Finally, 
Appendices [A] and [5] provide details of the derivations of 
Eqs. ([7]) and ( 13 ) respectively. 



II. EQUATIONS OF MOTION 

In this section, we present the derivation of the sec- 
ularized (Lindblad) form of the density matrix equa- 
tions of motion. Our approach initially follows the same 
steps as in Refs. [TH and [T|)] However, in order to pro- 
duce an approach that preserves positivity of the den- 
sity matrix throughout the range of validity of perturba- 
tion theor y, we must additionally secularize the resulting 
equations! 21 * 22 ! For completeness, we include in this pa- 
per also the steps which follow Refs. Q3] and [THl In this 
section we outline the main steps of the derivation, and 
present further details in appendix [A] 

We consider a single quantum dot, driven close to one 
of its discrete transition frequencies by a laser pulse with 
a time-dependent amplitude and frequency. For simplic- 
ity we assume that the excitation is circularly polarized, 
so that only one of the exciton spin states is relevant and 
we may neglect the formation of biexcitons. Thus we 
model the dot as a two-level system, which may be in 
the ground state, |0), or the one-exciton state, \X). The 



H dot = A(t)s z - n(t)s x , 



(1) 



in the rotating-wave approximation, and in a frame ro- 
tating at the instantaneous driving frequency u)(t). Here 
A(t) = Eq — ui{t) is the detuning between the transition 
energy Eq and the instantaneous driving frequency. O(i) 
is the time-dependent Rabi frequency, corresponding to 
the amplitude of the driving pulse, s = a/2, where the 
Pauli operator a z = \X)(X\ — |0)(0| describes the occu- 
pation of the exciton state, while a x = |A)(0| + |0)(A| = 
s_|_ + s_ describes the electric dipole moment of the tran- 
sition. We refer to s z = a z /2 as the inversion. 

We focus on the effects of acoustic phonons, which are 

the domi nant dephasing mechan ism in r ecent Rabi flop- 
pingSEH and ARp experiments j 4 | 5 | 9 | io | They couple to 

the dot via the deformation potential coupling, 



H c = s z ^2(g q b q + g*b\), 



(2) 



where q labels the phonon wavevectors, b q (b q ) is a 
phonon annihilation (creation) operator, and g q is the 
coupling constant. The phonon effects are controlled by 
the phonon spectral density, J(w) = J2 q \9q\ 2 ^( u ~ 00 q)- 
We take the modeP^for a GaAs/InGaAs quantum dot: 



J(w) 
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(3) 



with A = 11.2 fs K -1 , froJc = 1.4 meV. In Eq. (|| the 
low-frequency behavior J(ui) cx ui 3 arises from the cou- 
pling and density-of-states for acoustic phonons, while 
the high-frequency cut-off at w c arises from the size of 
the dot; confined excitons do not couple effectively to 
phonons of wavelengths smaller than the confinement. 

In the limit that A(t),f2(i) vary slowly with time, we 
may treat them effectively as constant and use the ap- 
proach discussed in Refs. [14] and [TBI so that the effect 
of the acoustic phonons can be found by transforming 
to the interaction picture, and using the Born-Markov 
approximation there. This requires that the interaction- 
picture density operator is approximately constant over 
the correlation time of the phonon bath (~ 1Aj c ), so that 
pit') — p{t) on the right-hand side of Eq. (A4). Equiva- 
lently, this means that the bath density of states should 
be flat over the effective linewidth of the system, as illus- 
trated in Fig. [I] Such an approximation is valid as long as 
neither the decay rate nor sweep rates (A/A, Cl/Vt) are 
too high, as both contribute to the effective linewidth. 
Note that applying the Born-Markov approximation di- 
rectly in the Schrodinger picture requires additionally 
that the density of states is flat on the scale of the range 
of position of the line, i.e., the energy scale of the Hamil- 
tonian, and this is not the case here. 

It is convenient to introduce rotated spin operators 
r = Rs, with R a rotation by angle tan -1 f2/A around 
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FIG. 1. (color online). Illustration of timescales required for 
validity of (time-dependent) Markov approximation, shown in 
the frequency domain. Solid red line shows the frequency de- 
pendent decay rate, illustrated by the phonon spectral func- 
tion J(uj). At low temperatures and slow driving (dashed 
line), the linewidth is sufficiently small that the decay rate 
does not vary significantly across the linewidth (bath correla- 
tion time is short compared to decay time). At higher tem- 
peratures, or faster driving, the linewidth grows (dot-dashed 
line) so that the decay rate does vary (decay time or sweep 
time become comparable to bath correlation time). The line- 
shapes illustrated are Lorentzians corresponding to transition 
rates T^ 1 =0.2 and 1.0 ps -1 , which may be compared with 
Fig-! 



the y-axis, so that the instantaneous system Hamilto- 
nian becomes H^ot — Ar 2 where A = \JVl 2 + A 2 . For 
the acoustic phonon coupling considered here, this yields 
(see Appendix |A| for further details): 



P = -PQ~p + QpP + PpQ^ - pQtP 
where P, Q are time-dependent operators of the form 

n 



2A 



(4) 



(5) 



Q(t) = / dvJ{y) \ dt'P{t') 



(6) 



and n v is the phonon occupation function at frequency 
v. After undoing the transformation to the interaction 
picture, this gives the density matrix equation form cor- 
responding to the results in Ref. [THl 

This equation is not of Lindblad form, and conse- 
quently it can lead to density matrix evolution which 
violates positivity. For the relatively short pulses in 
Ref. [T31 and [Tj)l one may readily check that this is not 
a problem. However, for our application to ARP pulses, 
positivity violation can occur at late times under con- 
ditions where the perturbative approximations required 
for Born-Markov remain valid; this is discussed further 
in Sec. 



This issue of positivity violation was discussed exten- 
sively in, e.g., Ref. [2TJ where it was shown that there 
exists more than one form of Markovian density ma- 
trix equation which faithfully represents the infinites- 
imal increment of the full density matrix evolution in 
the Markovian (perturbative) limit. However, although 
these different forms are equivalent regarding infinitesi- 
mal timesteps, those equations that are not of Lindblad 
form do not conserve positivity. A Lindblad form can 
nonetheless be derived by averaging so as to remove the 
rapidly oscillating terms in Eq. Q — such a procedure, 
known as secularization, yields a Lindblad form that in 
the perturbative limit is equally valid to Eq. Q. Further 
details are presented in the appendix; the result is: 

P = - [7o(f2) A)/2](r_r + p + pr_r + - 2r+pr_) 

- [ 7e (fi, A)/2] (r+r_p + pr + r_ - 2r_pr+) (7) 
-i[r z ,p](A + AE(Q,A)), 

where we have made explicit the time dependence of the 
decay rates due to the dependence on the slow variation 
of the parameters fi(t), A(t). This is a time-dependent 
generalization of the standard forrrP^ obtained from the 
secularized Born-Markov approximation in the in terac- 
tion picture, as used in some related contexts! 23 * 24 ! In ad- 
dition to now preserving positivity, it makes explicit the 
origin and nature of decay terms; the damping appears 
as a Lindblad form describing transitions between the in- 
stantaneous dressed states, with the phonon absorption 
and emission rates 
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9. 



(8) 



7e = 2 l2Ai ^ J ( A )K A ) + 1 ]- ( 9 ) 

These rates are shown in Fig. [2] for the spectral function 
Eq. pj. In addition, Eq. Q includes a phonon Lamb 
shift (Fig. |3| : the energy splitting of the dressed states 
is now A + AE with 



AE 



2A ) 



2A 



J{v) coth(i//2fcT) 
v 2 - A 2 



dp. 



(10) 



Figure [2] may be used to establish the validity of 
Eq. |7]) , by comparing the decay rates with the extent of 
their frequency dependence. We see that at the highest 
temperatures shown the peak damping rates, and there- 
fore the linewidths, can become a significant fraction of 
the width of the spectral function, as illustrated in Fig.[T] 
In this regime the Markovian approximation that p(t) 
varies slowly on the timescale 1/uj c breaks down, and a 
quantita tive analysis requires the solution of a non-local 
cquationPiSI However, as can be seen from the linewidths 
in Fig. [l] we expect qualitatively reasonable results over 
much of the parameter regimes shown, and the approx- 
imations should be quantitatively accurate in the low 
temperature regime, below 20K, typical of most coherent 
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FIG. 2. (color online). Phonon absorption (top) and emission 
(bottom) rates, Eqs. Q and Q, divided by the squared ratio 
of Rabi splitting to dressed-state energy splitting, Q 2 /A 2 , as a 
function of the dressed-state energy splitting A = \/Q. 2 + A 2 ; 
for resonant driving A = and the scaling is one. 
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FIG. 3. (color online). Phonon-induced change in the dressed- 
state splitting AE [Eq. ( 10 1] divided by the squared ra- 
tio of Rabi splitting to dressed-state energy splitting, i.e., 
AE/(Q 2 /A 2 ), as a function of the dressed-state energy split- 
ting A. Curves are temperatures IK (solid), 10K (dashed), 
50K (dot-dashed). 



control experiments. We note that the Markov approxi- 
mation amounts to approximating the spectral function 
with its constant value at the dressed frequency A. At 
high temperatures the damping at any one splitting will 
sample a finite range of the spectral function, so that 
we expect a weaker dependence of the effective linewidth 
on the energy splittings than indicated here, as well as 
non-Lorentzian emission lines. Similarly, we expect the 



Markov approximation to overestimate the Lamb shift 
and its frequency dependence at high temperatures. 

As well as becoming invalid at high temperatures, 
where the effective linewidth becomes large due to scat- 
tering, the approximations used above also fail if the time 
dependence of the parameters A(£),f2(i) becomes too 
strong. This is due to the finite bandwidth 1/T c hi rp ~ 
A/ A arising from the time dependence of the parame- 
ters. Alternatively one may understand this as arising 
from the fact that for small enough r c hi r p, the bath cor- 
relation time no longer is the shortest timescale in the 
problem. 



III. INVERSION IN ARP 

The above results apply in general to any time- 
dependent pulse sequence. We can in particular consider 
pulses corresponding to adiabatic rapid passage, in which 
the detuning A(£) is swept smoothly through zero with 
the intention of creating a one-exciton state. Under such 
a pulse there is an avoided crossing between the zero- 
and one- exciton states of Eq. ([I]) generated by the driv- 
ing field, and we aim to adiabatically follow this energy 
level, evolving from the initial ground state to the one ex- 
citon state. Acoustic phonon effects in this process have 
already been explicitly considered in Ref. [5] and [TUl for 
fixed-bandwidth pulses of the form 



A(t) = - 



at 



2(a 2 



OfA Qq / \ (11) 

W = / eXP ~ 1f 2 , 4\ 



where 0o is the area of the bandwidth-limited pulse be- 
fore the chirp is applied, and a is the spectral chirp 
The results obtained from Eq. Q for this form of pulse, 
plotted on a similar scale to those in Fig. 2 of Ref. H3 are 
shown in Fig. |4j As expected from the discussion above, 
the results are very similar at low temperatures and slow 
sweep rates, where the conditions for the Born-Markov 
approximation are well satisfied. 

As discussed in Ref. [9l the asymmetry about the line 
a = at low temperatures arises because absorption pro- 
cesses can be neglected j a <C 7e for T <C lj c , and as is 
clear from Eq. |7]) , emission can only occur when the dot 
is in the higher-energy dressed state. For a < the ARP 
protocol attempts to follow this higher energy dressed 
state, so that the process is susceptible to phonon emis- 
sion, whereas for a > it is not. The simplicity of Eq. |7]) 
allows one to further see that for the range plotted, the 
values of A(t = 0) = f2(i = 0) lie below the peak of scat- 
tering rates (Fig. [2]), hence the decrease of inversion with 
increasing pulse area visible within Fig. [4j For larger 
pulse areas, the central value of A can exceed this peak 
(at O w 6tt for a = 0), and inversion then increases with 
pulse area. 
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FIG. 4. (color online). Final exciton occupation probabil 
ity following an ARP pulse of the form given in Eq. (Ill 



calculated using the time-dependent Lindblad form given in 
Eq. j7|. Calculated using ro = 2 ps, at temperatures IK, 20K 
and 50K (top to bottom) for direct comparison with Fig. 2 of 
Ref. [9] For low temperatures, and slow sweeps, the results are 
very similar; for larger temperatures and higher chirp rates 
the limited Markov approximation used in deriving Eq. Q 
becomes invalid. 



A. Optimization of ARP 

The relatively lightweight effort of simulating Eq. |7]) 
allows one to rapidly investigate the effects of other po- 
tential pulse shapes and ARP protocols beyond that in 
Eq. (111. In the absence of decay, the question of how 
final 



et al.^. By considering the leading order non-adiabatic 
effects, they showed that these were minimized in the 
case where A(t) was independent of time. This implies 
that in the limit of large pulse areas (deep in the adi- 
abatic regime), maximum excitation should be reached 
when the chirp rate is adjusted to match this condition. 



The time dependence of A,Q given in Eq. (11) can 



not achieve a time independent A. Instead, other pulse 
shapes need to be considered, such as: 



A(i) = -A tanh 



f2(i) = f2o sech 



(12) 



which has a pulse area O when flo = Q/t. As discussed 
by Guerin et a/.^, the optimum condition f2 = A 
arises from the convergence of lines originating from the 
maxima of the Rabi oscillations. This is shown in the 
inset of Fig. [5] However, although such pulses are opti- 
mal in the isolated case, the differences in excitation near 
this line are exponentially small, and entirely dwarfed by 
the effects of phonon induced dephasing. In the presence 
of dephasing, not only are the sharp, exponentially small 
features washed out, but the optimum chirp rate a moves 
to significantly larger values, due to the reduction of the 
dephasing rate with at sufficiently large A. 




4 6 

Pulse area ©/7t 

FIG. 5. Final exciton occupation probability for the sech 
pulse [Eq. ( 12 1] with £7q = 0/r, Ao = olt for r = 6 ps at 



excited state population can be optimized for a 
given pulse area has been extensively studied by Guerin 



a temperature 4 K. The gray scale and contour labels show 
log 10 (l — P) where P is the final occupation probability. In- 
sets show the case without phonon induced dephasing, as in 
Ref. 1261 on the same gray scale (from -2 to 0, left inset), and 
on a gray scale with a larger range (from -10 to 0, right in- 
set). Main panel shows the effects of dynamical dephasing 
on the same pulse shape. Dashed line indicates the condition 
f2o = Aq which gives the optimal transfer in the absence of 
dephasing. 
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B. Thermalization enhanced inversion 

The discussion of the effects of acoustic phonons so far 
has been in terms of their reducing the final state inver- 
sion as compared to near-perfect inversion achieved deep 
in the adiabatic regime. There exists however a signifi- 
cant range of experimental conditions for which coupling 
to phonons can instead enhance the final state inver- 
sion. This effect has recently been discussed by Rciter 
et al.^H in the context of compensating for detuning of 
quantum dots. Even in the absence of detuning, coupling 
to phonons can enhance the final state inversion. 
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[shown in Fig.[6[a)] exists throughout the region between 
the solid and dashed lines in Fig. [6^c) . For large chirp 
rates and weak pulses, the inversion without coupling to 
phonons is already poor, and coupling to phonons in- 
creases the inversion, this corresponds to the behavior 
above the dashed line in Fig. [6^c). 

The origin of this enhancement at large q can be un- 
derstood by considering that for large q, the quantum dot 
state will come to thermal equilibrium with the phonon 
bath. Since the coupling to phonons depends on the pref- 
actor g(£l/A) 2 , the coupling to phonons will eventually 
switch off as — > 0. However, the larger the value of q, 
the later this switch off occurs, and so the longer system 
maintains thermal equilibrium with phonons. As the de- 
tuning A continues to increase at late times, the inversion 
of this equilibrated state therefore increases with increas- 
ing q. 

IV. PHONON SPECTROSCOPY 

A second application of the relative simplicity of 
Eq. Q is to see how the phonon density of states can 
be recovered from spectroscopy using an appropriately 
designed pulse sequence. This would in principle allow 
direct experimental confirmation of the model p honon 
coupling J(uj) oc u> 3 exp(— lu 2 /2u> 2 ) as widely usecPSEOD 
in modeling quantum dots. In the following we present 
and compare two approaches to this, based on either 
short-time or long-time behavior, incorporating sponta- 
neous decay in the long-time process. 



FIG. 6. (color online). Dependence of inversion on strength 
of phonon coupling in ARP. Panel (a) indicates how for a spe- 
cific value of a,Qo, there is a non-monotonic dependence of 
inversion upon phonon coupling. Panel (b) indicates the time 
evolution of s z at two values of q, relative to its ideal value 
for adiabatically following the dressed states. Panel (c) maps 
where non-monotonic dependence on q arises. This is deter- 
mined by two quantities: the grayscale indicates the value 
of Ma,xo<q<2{ds z /dq), and the solid red line is the bound- 
ary where this is zero. The blue dashed line indicates where 
ds z /dq\ q= o = 0. Between these lines, non-monotonic depen- 
dence as seen in panel (a) occurs. The magenta cross indi- 
cates the conditions used for panels (a) and (b). T = 4K, 
r = 5.68ps 

Figure [6] illustrates this potential enhancement by 
showing how the final inversion is affected if the decay 
rates are rescaled by a factor q, i.e. J(uj) — > qJ{uj), con- 
sidering a pulse shape: 

fl = Q sech (-) , A = -at. (13) 

For a wide range of parameters 0, a the dependence on 
q is non-monotonic: small coupling to phonons decreases 
the inversion, but further increase in coupling then in- 
creases final inversion. This non-monotonic behavior 



A. Modified ARP protocol 

The short-time approach uses a modified version of 
an ARP-like pulse sequence. We consider an ARP pro- 
tocol divided into pieces, so that the sweep of A(t) is 
interrupted by a wait period T w at a value A w before 
completing the ARP sweep. The initial and final sweeps 
serve to map between initial or final eigenstates of s z 
and eigenstates of r z during the wait period. The final 
inversion depends on the effect of the phonons during the 
wait time, sampling the phonon density of states at a fre- 
quency A w = y / A~f^+T£^. During this wait time, Eq. ^ 
reduces to pi = -j a pi + 7ePt.Pt = 7qPT ~ 7eP; in terms 
of the diagonal elements in the r basis. The deviation 
from inversion of the final state after a given wait time 
T w is thus given by: 

5-'.-^( 1 -«- h - +7 - )T -). <») 

where j aie = "f a ^ e (n w , A w ), and the numerator of the 
right hand side is j a for a forward sweep (A decreasing 
with time) and j e for a reverse sweep (A increasing with 
time). For long wait times, the excitations reach thermal 
equilibrium with the phonon bath, as expected^ and so 
s z becomes independent of the phonon density of states. 
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For short wait times (j a +^ e )T w <C 1, one finds | — s z ~ 
7o,e(^«jj 0)Tk,, thus by varying O^, one can directly map 
out the damping rate. 

In order to extract the phonon density of states with 
some accuracy, the pulse sequence must be carefully 
chosen. At low temperatures, 7a -C 7e , and so spec- 
troscopy using the forward sweep is hard to achieve — 
the reduction in inversion is tiny for times T w such that 
Twija + 7e) < 1 and tends to be dwarfed by effects of 
non-adiabaticity. Using a "reverse" ARP pulse produces 
a clearer signal. However, in order to have the signal 
dominated by the waiting time it is necessary for the wait 
time to be longer than the sweep, and the sweep to be 
sufficiently slow for all values of Q w . For this to be com- 
patible with T w -f e < 1 it is helpful to choose A w ^ so 
that the rate 7 e is suppressed by a factor (yi w /K w ) 2 < 1. 
Combining these considerations, the pulse sequence: 



A(i) = A w + A Q 



tanh 



t + t A 



tanh 



t-t A 



fi(t) 



tanh 



t + tg 

TO 



— tanh 



t-tg 



(15) 



gives the results shown in Fig. [7] with parameter values 
given in the figure caption. 




A during wait (ps)~ 



FIG. 7. (color online). Phonon spectroscopy using divided 
ARP pulse. Main figure, deviation from inversion following 
the protocol in Eq. flXSj) (solid), compared with the approx- 
\ (das 



imation of Eq. |14j) (dashed line), plotted by varying Q w . 
Inset, time dependence of A(t) (solid) and Q(t) (dashed), 
with parameters A w = 2ps , An = — 10ps _1 ,fA = T w /2 = 
3ps,rA = 1.2ps,tn = 15ps, tq = 5ps, T = 4K. 

There is a reasonable match between the prediction of 



Eq. ( 14 ) for decay during the wait time and the actual 



inversion, but the match is not perfect due to the effects 
of phonons during the initial and final sweep, as well 
as some remaining degree of non-adiabaticity (at small 
flu,). Nonetheless, one may invert Eq. (14 1 in order to 



extract the phonon density of states from the measured 
inversion, and the result of this procedure is shown as 
the dot-dashed line in the inset of Fig. [8] The reason- 
able agreement confirms that for such a pulse the final 



inversion contains sufficient information to extract the 
phonon density-of-states; in practice one might include 
corrections to Eq. ( 14 ) by comparing an experiment di- 
rectly to Eq. 0. 



B. Steady state of driven open system 

An alternate approach to reconstructing the phonon 
density of states arises by considering the long time be- 
havior, allowing for both spontaneous decay as well as 
coupling to acoustic phonons, i.e. 



P - ( K / 2 )( s +S-P + ps+s- - 2s-ps+). 



(16) 



In particular, we consider the long time behavior under 
a constant driving, f2(t) = f2o, A(t) = Ao- In this case, 
the long time behavior should be understood as a steady 
state arising from the balance of coupling to the phonon 
reservoir and the decay due to coupling to a photon reser- 
voir. The pumping term Q is necessary only in order to 
enable the phonon coupling to create and destroy excita- 
tions of the quantum dot. In the absence of spontaneous 
decay terms (no photon reservoir) , the steady state would 
trivially be the thermal equilibrium of the system Hamil- 
tonian, due to equilibration with the phonon reservoir. 
In this case, the phonon density of states does not ap- 
pear, only the phonon temperature matters. However, 
the coupling to the photon reservoir drives the system 
into a non-thermal steady state, in which the balance of 
spontaneous decay rate and phonon coupling determines 
the state. One may then read out the phonon density of 
states from the steady state inversion achieved. 

In the limit of small decays, i.e., n,j a ,j e <C Q,A, 
one may analytically solve the steady state equations for 
s x,y,z by expanding the steady state equations with re- 
spect to the decay rates (see appendix [B] for further de- 
tails). One finds that at leading order: 



— A 2 k + 2AA( 7n - 7e ) 
(ft 2 + 2A 2 )K + 4A 2 ( 7a +7 e )- 



(17) 



As anticipated, this expression involves the ratio of 
phonon and spontaneous decay terms. In the limit 
k — > 0, one recovers the thermal equilibrium result so 
that s z = (-A/2A)( 7e - J a )/{7a + 7e), in which case 
the phonon density of states cancels, and no information 
about the phonon bath (other than temperature) appears 
in s z . However for k ^ 0, the steady state depends on the 
ratio of 7° / ft,7 e /k, and this in turn allows the phonon 
density of states to be extracted from the final state. 

Figure |H] shows how the phonon density of states can be 
reconstructed by extracting the steady state for a fixed 
value of f2, and varying the pump detuning A. By in- 
verting Eq. (17), and assuming the phonon temperature 



is known, one may extract the effective phonon density of 
states as shown in the inset. The density of states recon- 
structed this way matches the actual density of states 
used in the density matrix evolution very closely. It 
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FIG. 8. (color online). Main figure: Steady state value of s z 
with driving at Q, — O.l(ps) -1 , with a spontaneous decay rate 
of k = 2(ns) _1 as a function of the detuning, at T — 20K. 
Inset: comparison of actual phonon density of states J(u>) 
with the values reconstructed by inverting E q. (|17[ ) and from 
the modified ARP approach, inverting Eq. ( |14| (parameters 
for modified ARP as in Fig. ffh. 



may also be possible to use recently developed varia- 
tional approaches^ to extend such phonon spectroscopy 
to more strongly coupled systems. 



V. LINDBLAD VS NON-LINDBLAD 
APPROXIMATIONS 

As mentioned above, the question of which Marko- 
vian approximation best corresponds to a given full den- 
sity matrix equation was discussed extensively by Dum- 
cke and Spohn^D. The conclusion there is that in the 
limit of short bath correlation times, there exist multi- 
ple Markovian approximations which give have the same 
order of validity, as defined by the limiting behavior of 
U^approx _ pfuii|| (where || . . . || is the trace norm) as cou- 
pling to the bath vanishes. These different Markovian 
approximations are related by whether or not one aver- 
ages over fast oscillating terms in the interaction picture. 
Without such averaging we reach Eq. Q, while averag- 
ing leads instead to Eq. ([7|. However, positivity of the 
density matrix is only preserved for the Lindblad form in 
Eq. (m), and so Dumcke and Spohn^ conclude that only 
this approach is correct. 

It is often assumed that such a point is irrelevant as 
small decay rates (as are required for validity of the 
Markov approximation) imply that any possible viola- 
tion of positivity is negligible. However, for problems 
involving weak decay and long time evolution, such as 
the current problem, positivity violation can occur for 
Eq. @ even in regimes where the Markov approximation 
appears to be valid. This leads to unphysical results, as 
shown for example in Fig. [9j We note that for the Rabi 
oscillations studied by Ramsay et alJ^ no positivity vi- 
olation occurs for the parameters used and the results of 



Eq. Q and Eq. (JtJ) are hardly distinguishable. Nonethe- 
less, Eq. Q is required to produce valid results across 
the problems considered in the current paper. Further- 
more, as shown in the inset to Fig. [9j the Lindblad form 
continues to gives qualitatively plausible results even at 
higher temperatures where Eq. Q very strongly violates 
positivity. 
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FIG. 9. (color online). Comparison of the dynamics obtained 
from the Lindblad [Eq. (JtJ , dashed curves] and non-Lindblad 
[Eq. Q, solid curves] forms of time-dependent Markovian ap- 
proximations, for the Gaussian ARP pulse, Eq. (Ill, with 
to = 2 ps, T = 1 K, a = 30 ps 2 , 0o = 5ir. For each ap- 
proach both the magnitude of the psuedospin (Bloch) vector 
|s| (left axis) and inversion s z (right axis) are shown. At this 
low temperature, as discussed above, the Markov approxima- 
tion holds reasonably well, yet the non-Lindblad form leads 
to unphysical results jsj, \s z \ > 0.5. The inset shows the in- 
version for a pulse with a — 36 ps 2 , Qo = tt at T — 65 K; 
here the non-Lindblad form rapidly gives unphysical results 
\s g \ > 0.5, while the Lindblad form gives qualitatively reason- 
able dynamics. 



There have been suggestions that time-dependent 
density-matrix equations which are not of the Lindblad 
form, such as those derived from the time-convolutionless 
(TCL) approach,^ may in some cases ensure positivity!^ 
However, Eq. Q is essentially the result of the time- 
convolutionless approach at second order (TCL2). The 
only difference is that for TCL2, the lower limit on in- 
tegrals over t' such as Eq. ^ is t' = 0, rather than 
t' = — oo, corresponding to starting the system at t = 
in a factorized state. Even for time independent A, f2 
this means the coefficients in the density matrix evolution 
become time dependent, but eventually decay toward a 
steady state value. As such, for the ARP pulse, as long 
as the pulse duration is long compared to the bath corre- 
lation time, these additional time dependencies die out, 
and TCL2 becomes equivalent to the positivity-violating 
form in Eq. Q. 
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VI. CONCLUSIONS 

In conclusion, we have shown how the Bloch-Redficld- 
Wangsness theory may be used to derive a time- 
dependent Lindblad form describing the dephasing of 
quantum-dot excitons by acoustic phonons in the pres- 
ence of a driving laser field. We have outlined the ap- 
plication of this theory to recent ARP experimental on 
single quantum dots and predict, in agreement with nu- 
merical work,® that phonons have a pronounced effect on 
ARP even at cryogenic temperatures. Their effect can, 
however, be almost eliminated by an appropriate choice 
of pulse shape. Furthermore, this pulse-shape depen- 
dence could allow forms of phonon spectroscopy based 
on ARP pulses or on off-resonant continuous- wave excita- 
tion. More generally, our approach captures the physics 
of dynamical, excitation-controlled dephasing, in which 
the driving field changes the eigenspectrum of the dot, 
and hence the decoherence and scattering rates. The 
resulting Lindblad form is straightforward to simulate, 
gives qualitatively reasonable results over wide parame- 
ter regimes, and is expected to be quantitatively accurate 
at low temperatures for slow pulses. It can be applied to 
a wide variety of pulse sequences, and the approach gen- 
eralized to address a wide range of problems relating to 
the decoherence of solid-state qubits, such as the lim- 
itations on the creation of entangled states in coupled 
quantum dotsJ^SHSl the persistence of entanglement in 
interacting solid-state systems, and the emission s pectr a 
of solid-state qubits in the strong-coupling regimePSESI 
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Appendix A: Derivation of secularized density 
matrix equation 

This appendix provides further details of the deriva- 
tion of the time-dependent Markovian Lindblad form in 
Eq. Q. Following the usual approach pi working in the 
interaction picture, the effects of the system-bath cou- 
pling on the system density matrix can be included to 
second order by writing: 



dt'Tr, 



H c (t), H c (t'),p(t')®p B (Q) (Al) 



where tilde indicates the interaction picture, ps is the 
phonon bath matrix, and H c is given in Eq. (2}. 



In the case that A, fl vary slowly with time, we may 
effectively transform the coupling Hamiltonian to the in- 
teraction picture by using the instantaneous eigenstates 
(dressed states) giving 



H c (t) = P(t)$(t) 



el tH dotSze -UH dot ^2(g qbqe -^t + g* q b\e 1 ^), (A2) 



where v q is the phonon frequency. We transform to the 
instantaneous dressed states, defining spin operators r = 
Rs, with R a rotation by angle tan" 1 f2/A around the 
y-axis. Then the coupling operators become 



P(t) 



A 

A 5 



ft 
2A 



(r+e jA * + r_e- 4At ), (A3) 



where A = Vft 2 + A 2 is the splitting of the instanta- 
neous eigenstates, and r± cause transitions between these 
states. 

Note that for a time-dependent Hamiltonian, the true 
interaction-picture form is obtained with a unitary trans- 
formation involving a time-ordered exponential, U — 
Te~ l $ H ^t(t )dt ^ anc j ^his operator is not generally ap- 
proximated by the form in Eq. (A2). However, in the 



following we use the form in Eq. (A2) only to calculate 
the effects of phonons on short timescales, t c ~ l/w C! 
with the final equation for the density matrix obtained 
by undoing this formal transformation. Thus we expect 
that the Hamiltonian part of the dynamics is not approx- 
imated in the result, while the dissipative part is correct 
provided t c dA/dt,t c d£l/dt << A. These conditions are 
well satisfied for the protocols considered in this paper. 

With these forms of P(t), $(t), we may then follow the 
normal steps of tracing over the phonon bath to give the 
system density matrix equation: 



dvJ(v) J dt'{ [P(t)P(t')p(t') P(t')p(t')P(t)} 

x \(n v + l)e- iu{ t-^ + r^e^*" 4 ' 5 

-[P(t)p(t')P(t')-p(t')P(t')P(t)} 

(n„ + l)^*-*') + rve"^*-*')] }, 

(A4) 



where n v is the thermal occupation of the phonons at fre- 
quency v. After performing the integrals over frequency 
the remaining integral contains factors which are sharply 
peaked functions of t — t', decaying over a timescale 1/uj c . 
p(t') and A(t'), f2(i') vary little over this timescale, and so 
may be approximated by their values at t. However, P(t') 
may vary over this timescale due to the time-dependence 
arising from the transformation to the interaction pic- 
ture. If we approximate p(t') ~ p(t) and perform no 
other steps, this leads to Eq. Q, which as noted before, 
is not of Lindblad form. 

Following Dumcke and Spotm^, the corresponding 
Lindblad form arises by "secularizing" Eq. (A4|. This 



10 



corresponds to averaging the above equation over a time 
short compared to decay rates, but long compared to the 
timescales of the system Hamiltonian — the fact that 
such a timescale exists is implicit in the use of a per- 
turbative (Born) approximation. We start from Eq. Q, 
with P(t) as defined in Eq. and writing Qit) in the 
generic form: 



Q(t) = r z r z + r H 



iAI 



(A5) 



thus becomes: 
P= ^(r z + r;) (r 2 pr z - p) 



ft 

2A 

ft 

2A 



[r+ {r-r+p - r +/ 5r_) + r* (pr_r+ - r+pr_)] 
[r_ (r + r_/5 — r_pr + ) + T*_ (pr + r- — r_/5r + )] 



The vanishing of the phonon density of states as uj — > 
ensures that T z +T* Z — 0, and the remaining terms take 
the form of the Lindblad decay and phonon Lamb shift 
terms as given in Eq. 



Appendix B: Steady state of driven open system for 
small decay rates 



This follows directly from performing the integrals in 
Eq. ([6]), and so T Z ,T± are various frequency integrals 
over J(u). Multiplying P(t),Q(t) and integrating over a 
time long compared to 1/A, only those terms with equal 
and opposite t dependence will survive, i.e. those terms 
involving r z r z ,r + r_ or r_r + . The secularized Eq. Q 



This appendix gives further details of the derivation of 
Eq. (17), the steady state s z found for the driven system 



in the presence of both phonon induced dephasing and 
spontaneous decay (photon emission). Including spon- 
taneous decay as in Eq. ( [16] ) in addition to the density 
matrix evolution given in Eq. ([7| one may write equations 
of motion for the Bloch vector s x y g of the form: 



(7a - 7e) 



2ft' 



A 2 



(7a + 7e) 



Sy 



(«/2 + 7a + le)Sy + 



k A , Aft, , 

2 A" le ' + A2"^ 7a + le > S * 



fts,, 



AS V + ~j2^ a +le)S z 

2A 2 + ft 2 , 

(7a +7e) 



A 2 



(Bl) 
(B2) 
(B3) 



In the limit of vanishing decay terms K,7 a ,7 e — > 0, it s x = — (ft/A)s z as the solution of s v = up to order 

is clear that the steady state requires As x + fts z = 0(k 2 , 7 2 , 7 2 ). Using this, the remaining two equations 

0, s y = 0. At leading order in the decay terms, one s x>z = can be solved by eliminating s y , to give Eq. ( 17 1, 
finds s y ~ 0(ft,7 a ,7 e ), and so one may continue to use 
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